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We apply the Lanczos method to a 2-leg ladder with mixed spins of magnitudes (Si, S2) = (1, 1/2) 
located at alternating positions along the ladder. The effect of dimerization 7 is also considered 
according to two different patterns. A Spin Wave Theory (SWT) is applied to this model predicting 
one gapless branch with ferromagnetic properties and another gapful branch with antiferromagnetic 
nature as low energy excitations of the model. We compute the ground state energies, Ferro- and 
AF-excitation gaps, magnetizations and correlation functions as a function of J' and 7 which results 
into a fine estimate of the phase diagram. The Lanczos results are compared with the SWT analysis 
and a qualitative agreement is found but with numerical discrepancies. We also study numerically 
the Spin-Peierls instability and find that it is conditional for any value of J' £ (0, 2) and both 
dimerization patterns. 

PACS number: 76.50.+g, 75.50.Gg, 75.10.Jm 



I. INTRODUCTION 

Mixed spin chains with alternating spins (1,1/2) are 
paradigmatic examples of strongly correlated systems 
exhibiting excitations of different types, namely, both 
gapless and gapped excitations in the low energy spec- 
trum |l] ^] . This is a manifestation of their ferrimag- 
netic character and they are responsible for their unusual 
properties. On the contrary, one-dimensional integer- 
spin Anti-ferromagnetic Heisenberg (AFH) model has 
a unique disordered ground state with a finite excita- 
tion gap while the half-integer spin AFH chain is gap- 
less [pi. These ferrimagnetic quantum chains have re- 
cently attracted much attention and several techniques 
have been devoted to their study : analytically, such as 
Spin Wave Theory (SWT) Pjj| , Quantum Renormaliza- 
tion Group H, Variational [|l|, and numerically, such as 
conventional |§, transfer matrix DMRG @, QMC |] 
etc... Moreover, the interest is increased due to the fact 
that experimental realizations of these one-dimensional 
ferrimagnets also exist. For instance, the oxamato com- 
pound NiCu(pba)(D20)3 ■ 2D 2 exhibits a structure of 
a (jSijSa) = (1, 1/2) ferrimagnetic chain j?J where the 
spins 1 are located at the Ni sites and the spins 1/2 at 
the Cu sites. It is found that below T = IK it exhibits 
Long Range Order (LRO). Materials realizing other mag- 
nitudes of spins (Si, S2) have been also synthesized [||. 
In all these compounds it is assumed that the interchain 
coupling is much smaller than the intra-chain interactions 
so that the system becomes a set of effectively uncoupled 
quantum chains. 

On the other hand, many interesting investigations 
have been devoted to spin ladders. They consist of cou- 
pled one-dimensional chains and may be regarded as 
interpolating one- and two-dimensional systems. The- 
oretical studies have suggested two different universal- 



ity class for the uniform spin ladders, i.e. the anti- 
ferromagnetic spin-i ladders are gapful or gapless, de- 
pending on whether ni (the number of legs) is even or 
odd [pi. Some experimental observations are in agree- 
ment with this results |h]] . However mixed-spin ladders 
may change this universality. In a recent work p"T[ we 
have considered a ferrimagnetic 2-leg ladder and conjec- 
tured its phase diagram by using an approximate Quan- 
tum RG analysis. We have shown that for the homoge- 
neous antiferromagnetic couplings the model is gapless 
but an energy gap may appear in a specific type of con- 
figuration of the coupling constants. 
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FIG. 1. A picture of the dimerization patterns considered 
in a 2-leg mixed spin ladder: a) Columnar Bond Alternation 
(CBA), b) Staggered Bond Alternation (SBA). Small circles 
are S — \ and large ones are S = 1. 

In this work we present the results of an extensive nu- 
merical analysis using the Lanczos method [ fl0| applied 
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to a large class of ferrimagnetic 2-leg ladders for the first 
time (see f| for a review on ladder materials.) Apart 
from accurate numerical results which show the effect 
of inter-chain interaction in a ferrimagnetic system, we 
obtain the phase diagram of Ref. [jll| with more accu- 
racy. In the present phase diagram, which is obtained by 
Lanczos results, we specify the phase boundary between 
the gapless and gapful phases more accurately and show 
that the transition between two gapless phases appears 
at zero inter-chain coupling and is first-order. An an- 
alytic Spin Wave Theory is also introduced to describe 
the low-energy excitations of this model. Moreover we 
will consider the effect of bond-alternation on the Spin- 
Peierls instability and its conditional or unconditional 
character. The model we have considered here, depicted 
in Fig.l, can be assumed as the first step to consider the 
inter-chain interactions in a ferrimagnetic material. (For 
instance, MnCu(pbaOH)(H20)z has such kind of struc- 
ture |T^| .) In our model mixed spins have lengths S\ = 1 
and S% = 1/2 and are located at alternating positions in 
this bipartite lattice, making up the sub-lattices A and B, 
respectively. One may think of two (1, 1/2) ferrimagnetic 
chains which interact via rungs in a ladder realization 
where different kind of spins are located on each rung 
(see Fig.|l|). Specifically, the isotropic Heisenberg Hamil- 
tonian we consider is the following ( with i — 2n + a) , 

2 --i 

ff = J EE + S^-S^ + l) 

+ [1 + 7 (Q) S^ a) {i + l)-S[ a \i + 2)) 

+ J'J2 [ S i 1} ( 2n + 1 )- s 2 2) (2n + 1) 

71 = 

+ S^ 1) (2n + 2)-Sf ) (2n + 2)] , (1) 

where S^(n) denotes the quantum spin-1 at site n in 
the leg a = 1,2 of the ladder, similarly S 2 a \n) for 
the spin-1/2, J and J' are coupling constants along 
the legs and the rungs respectively, and 7 < > Q - ) (n) is the 
dimerization patterns. We consider two different dimcr- 
ization patterns : i) CBA (Columnar Bond Alterna- 
tion, FigJ^(a)) for which ^ a \n) = (-l) n+1 7 0, and 
ii) SBA (Staggered Bond Alternation, Fig.|l|(b)) where 
7 ( Q )(n) = (-l) Q + 1 +" 7 ||@. In order to keep the sys- 
tem always in the antiferromagnetic regime, we restrict 
the range of variation as \ j\ < 1. We use periodic bound- 
ary conditions along the legs of the ladder. The number 
of sites is N = 2 x L, where L is the length. 

The paper is organized as follows: in next section 
we introduce the spin wave theory applied to a non- 
dimerised (7 = 0) ferrimagnetic (Si,^) ladder. Dis- 
persion relations, ground state energy , energy gap and 
magnetization of this model are obtained within SWT in 
order to explain the nature of elementary excitations in 



this system. In Sec. Ill we present our numerical results 
obtained with the Lanczos method and we contrast them 
against the previous SWT results. In Sec. IV we extend 
our Lanczos results to include a numerical analysis of the 
phase diagram for a 2-leg mixed spin ladder to specify the 
boundary between the gapless and gapful phases. In Sec. 
V. the spin-Peierls instability is discussed by computing 
the magnetic energy gain of the present model. We fi- 
nally present our conclusions in Sec. VI. 

II. SPIN WAVE THEORY 

In this section we shall perform a mean-field treatment 
of the ferrimagnetic ladders and for simplicity we shall 
concentrate on the case of zero dimerization 7 = 0. 

The application of the Lieb-Mattis theorem to this 
ferrimagnetic ladder predicts the existence of a Ground 
State (GS) with Total Spin S G given by S G = L(S 1 ~S 2 ). 
This leads to the existence of a non-vanishing value of 
the magnetization at zero temperature. Then, a natural 
issue arises, namely, the study of classical ferrimagnetic 
order v.s. quantum antiferromagnetic fluctuations ||. 
This system presents gapless Ferromagnetic excitations 
(Goldstone modes due to broken symmetry) with spin 
Sq — 1 and gapped Antiferromagnetic excitations with 
spin Sg + 1. 

Thus, as we expect an ordered ground state, it is natural 
to perform a spin wave theory analysis to describe the 
low energy quantum fluctuations around a classical ferri- 
magnetic ground state. We shall employ the SWT in the 
linear approximation. 

To this end, we first make a standard Holstein- 
Primakoff transformation on the (1, l/2)-spins of eq. (Q). 
Let be the modes in momentum space associated to 

the sub-lattice A for the legs a = 1, 2, and similarly 
for the modes in the sub-lattice B. The Hamiltonian can 
be written in the following form. 

H S WT = H rung + ^ H leg 
a=l,2 

H rung = —LJS1S2 

+ j'E^(4 (1) 4 15 + 4 (2) 4 2) ) + + bfh^)] 
+ j'v^£(« 2) + W + 4 (1) e + W) 

k 

< = -LJ Sl S 2 + 2jJ2(Sial^a^ + S^h^) 

k 

+ jv^E 2Cos ( fc )(4 a) ^ Q) +4 a) b[ {a) ) (2) 
A,- 

As it is apparent, this is not enough to diagonalize the 
Hamiltonian with a Bogoliubov transformation because 
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the degrees of freedom of each leg appear coupled to- 
gether. Fortunately enough, we can devise a second 
transformation by introducing a couple of symmetric 
s', s" and antisymmetric a', a" fields as follows: 



where , -B^ are the quasi-particle creation oper- 

ators and the frequencies w^'Qe), LOg(k) of the several 
excitations take the following expressions as a function 
of the model coupling constants, 



« (a) = ^(4 + (-i) 1+ v fc ), 
6 w = -L(4 + (_i)i+X) ! 



(3) 



J A 
,(± 



(fc) = i(d -C 2 + ^/(C 1 + C 2 r-4C 2 3 i± \k)) 
\k) = ~(C 2 -Ci + V(Ci+C 2 )2-4C 3 2(±) (fc)) 



with a = 1,2 represents leg index. Then, it is possible 
to show that the resulting SWT Hamiltonian takes the 
following form in terms of these fields: 

Hswt = -LS 1 S 2 (2J + J') + h(s) + h(a) , (4) 

where 

h(8) = C 1 y EsUk + C a y E*k% (5) 



+ V^E( J2cosfc + J ')Wl + *k*k) > 

k 

h{a) = C 1 Y,4< + C 2 Y,H^l (6) 

k k 

+ y/thJhJ2( J2cosk - J ')( a 'k a 'k + 4 a 'k) > 

k 

and now using a set of Bogoliubov transformations we 
will arrive at the following expressions for these Hamil- 
tonians h{s), h(a) as, 

h( S ) = ]r[c (+) w +4 +) w4 (+) 4 +) + - , 
A,- 

(7) 



C (±) (fc) = -^(Cr - C 2 ) + ^(C 1+ C7 2 )2-4C 3 2(±) (fc) , 
d = (2 J + J')Si C 2 = (2 J + J')S 2 , 

Cf\k) = \J SiS 2 (2J cosk ± J') , (8) 

In Fig.|] we have plotted the four branches of the dis- 
persion relations (d^ (k),Wg(k). For vanishing inter- 
chain coupling J' = (Fig. ^ (a)) we reproduce the 
ld-chain result. As the interaction is switched on, we 
appreciate the splitting of each zero-coupling branch into 
another two new branches in the case of ladder. The ap- 
pearance of four modes is related to defining two boson 
operators for each leg of ladder. At fc = we define the 
gap as gi = \ojg \<S) — (0)|, which goes to zero with 
vanishing J' as can be seen in Fig.||. We also define the 
gap g 2 = (0) — (0)|, which remains finite VJ' 
(see Fig |H). So we consider g 2 as the actual energy gap 
which is defined by SWT. The lowest energy dispersion 
relation is J B +} which shows gapless excitations at k = 0, 



,(+) 



(fc -> 0) = 



2S\S 2 J 2 
|Si-S 2 | 



(9) 



The first excited branch which shows an actual gap for 
any value of J' is uj^ with the following behaviour in 
the long wave length limit, 




,(+) 



(fc^0) = (2J + J')\S 1 -S 2 \ 



2S\S 2 J 2 
|5i-5 2 | 



(10) 
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FIG. 2. Dispersion relations for a ferrimagnetic ladder for 
different values of rung coupling namely J' = 0.0, 0.2, 1.0, 2.0. 
The J' = case (a) is the same as one-dimensional ferrimag- 
nets. 



Eq.(||) is similar to a ferromagnetic dispersion relation, 
then we expect a ferromagnetic behaviour for these ex- 
citations. Moreover it can be shown that the oj^\k) 
modes are created by acting on the classical ferrimagnetic 
ordered state (|0» with S~ = £ n E Q (^ W (n) + 

S 2 ^ a \n))e lkn , where k 7^ 0. The outcomes of this oper- 
ation are those states which have S* ot — L\Si ~ S 2 \ — 1. 
These states are annihilated by S^ ot . Therefore we con- 
clude that these states have Stot = L\Si — S 2 \ — 1 which 
is another sign for the ferromagnetic behaviour of u>^ 
modes. The magnetization of Si = 1 sub-lattice shown 
in Fig.(^j) verifies the reduction of magnetization from 
its classical value ((Sf)ciass. = 1) which shows that the 
quantum fluctuations of the lowest branch has a ferro- 
magnetic property. Similarly, the oj^ modes are ere- 
ated by acting on |0) with S+ = ^£„ £ Q (S+ (Q) (n) + 
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S2^ a \n))e tkn , where k ^ 0. It is also possible to show 
that those states have Stot = L\Si — Sa| + 1 which is the 
property of an anti-ferromagnetic excitations, i.e. the ex- 
cited states have larger spin than the ground state. All of 
these explanations for the excited states are examined for 
a finite system size (N = 2L = 8, 12, 16, 20) by Lanczos 
method where we will discuss them in the next section. 

Using the quantities (^(k) and eqs. we have 

also computed the GS energy per site E /JN as follows, 

E Q = -NS 1 S 2 (2J + J') + ^(C (+) (fc) + C (_) (*))- (H) 

fe 

It is plotted as a function of J' in Fig|j. 
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FIG. 3. Lanczos v.s. SWT results for GS energy, as a func- 
tion of the interchain coupling J'. Inset similarly shows the 
Lanczos results v.s. SWT for energy gaps. 

In a similar way the sub-lattice magnetization m\ as a 
function of the couplings is calculated, 



m 1 = ^{Sf\i) + Sf\i)). 



(12) 



The other sub-lattice magnetization mi — (Sf ) can also 
be obtained by using the identity : m\ + m-i — Si — S 2 - 
We have plotted mi as a function of J' in Fig.[|. We will 
discuss on this result in the next section. 



III. NUMERICAL ANALYSIS: LANCZOS 
METHOD 

Let us now present our numerical Lanczos results. 
We have studied ladders with a number of sites N — 
2 x L with L = 4,6, 8, 10 due to the constraint of pe- 
riodic BC's. We always set J = 1 and vary J' = 
0.0, 0.2, 0.8, 1.0, 1.2, 1.8, 2.0. The dimerization parameter 



is varied from 7 = 0.0 to 7 = 1.0 in steps of 0.2 for both 
CBA and SB A patterns. 

As a testground, we have checked several quantities at 
J' = 0, like the GS energy per site eo and Antiferromag- 
netic gap Ao and we find excellent agreement with the 
DMRG results for the spin chain as shown in Table 1, 
namely, e (DMRG) = -0.727 and A„(DMRG) = 1.759 
g@. 

In Fig. we plot the variation with J' of the GS energy 
and we find that the extrapolated numerical results agree 
very well with the SWT results in the whole range. How- 
ever, for the Antiferromagnetic gap (see inset of Fig||) 
we find clear numerical differences with SWT. For com- 
parison, we have plotted the two gaps, 31,32 introduced 
before, because for the strong coupling region around 
J' ~ 2, both gaps turn out to be very close to each 
other. We clearly find that the exact gap obtained with 
Lanczos is well above the values obtained with the SWT 
approximation. 

Magnetization (m,) 




O © Lanczos 

A ASWT 



FIG. 4. Sub-lattice magnetization of spin Si = 1 ({Si)), 
Lanczos v.s. SWT results as a function of J'. 

Likewise, we have computed several sub-lattice mag- 
netization and in Fig. |^ we show mi = (Sf(l)) for the 
spin-1 case at site n = 1 and we have checked that it is 
the same at all other positions due to translational in- 
variance. We find that the SWT calculation essentially 
underestimates the Lanczos result. The discrepancies 
between SWT results for energy gap (inset of Fig. ||) 
and magnetization (Fig. |J) means that the linear SWT 
does not provide us with a complete account of quan- 
tum fluctuations around the classical ferrimagnetic or- 
der. However the ground state energy obtained by both 
methods are in good agreement. For instance in the ho- 
mogeneous case (7 = 0, J = 1, J' = 1) the SWT result is 
e (SWT) = ffr = -0.96475 which has only 0.7 percent 
error with respect to Lanczos result (Table. 1). 
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FIG. 5. Plot of the GS energy (Lanczos) v.s. 1/N for CBA 
and SBA dimerization patterns. Convergence to thermody- 
namic limit is very sharp for all dimerization value (7). 

We have also studied the effect of dimerization in the 
ferrimagnetic 2-leg ladder in the Columnar (CBA) and 
staggered (SBA) configurations. In Fig. |] we present the 
results of Lanczos calculations for J' = 1 and varying the 
staggering parameter 7 in the whole range. We first check 
that our results are well converged for the sizes considered 
(N=8,12,16,20). The ground state energy converges to 
the value of thermodynamic limit as rapid as in the one- 
dimensional ferrimagnets. It means that the finite size 
effects are small in the ferrimagnetic ladders as well as in 
1-D models. We then observe that the effect of dimeriza- 
tion (7) is generically to decrease the GS energy and that 
the columnar configuration produces a bigger energy de- 
crease than the staggered pattern. This shows that the 
CBA pattern is more stable than the SBA configuration. 
The stability is improved by increasing dimerization. 



for different values of J'(0.8, 1.0, 1.2, 2.0) which shows the 
same behaviour. We have also checked for some more val- 
ues of J' and the same results are obtained. Moreover, we 
have found that this also holds true when dimerization is 
present and in Fig.^|(c) we show this for the SBA configu- 
ration which shows that we have gapless excitations with 
ferromagnetic nature for all values of J' and the dimer- 
ization parameters (7). In Figs.^|(b),(d) we have plotted 
the AF-gap E q (Sg + 1) - E q (Sg) for some values of J' 
and 7. We find that the AF-gap remains finite in the 
thermodynamic limit, moreover the interchain coupling 
J' increases the gap both with or without dimerization. 
The notion of Ferro-gap and AF-gap is related to the ex- 
citations with spin lower or higher than the spin of the 
ground state, respectively. This is in complete agreement 
with the SWT explanation of the low energy spectrum. 
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FIG. 6. Plot of energy gaps (Lanczos): (a)-(b) Ferro-gap 
and AF-gap without dimerization v.s. 1/N for several cou- 
plings J', (c)-(d) Same with dimerization (SBA). 
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Table 1. Extrapolated (N — > 00) Lanczos results for GS 
density energy eo = Eq/JN and AF-gap Aq. 



The Lanczos analysis of the ferromagnetic and anti- 
ferromagnetic gaps are shown in Figsj^. In Fig.|^(a) we 
clearly find that in the absence of dimerization 7 = 0, 
the Ferro-gap E (Sg — 1) — E (Sg) scales to zero as the 
system size N goes to 00. This corresponds to the gapless 
excitations in the SWT analysis. We have plotted data 



We have presented the extrapolated values of ground 
state energy per site (eo) and AF-gap (A ) in Table 1. 
The scaling form which we have considered for the ex- 
trapolation to N — ► 00 is a power-law fit. We have im- 
plemented the following function for the scaling form of 
the ground state energy, 

e (N) = e (oo) + ^ (13) 

where a and v are constants and obtained to have least- 
square root error. We have examined a polynomial form 
but the best fitting was obtained by the power-law fit 
of Eq.( [l3|). In the case of energy gap we considered 
both Eq.(|l3|) and an exponential form, (i.e. Ao(A^) = 
A (oo) +B exp{— \c\N)/(N^)) where B, c and fi are con- 
stants. Both scaling forms give us the same results for 
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Ao(oo) and fit very well to the data of Fig.(||). Although 
from the SWT calculations we obtained the k 2 depen- 
dence for both gapless and gapped spectrum (Eqs.(||, |l~C|)) 
which impose that the scaling form should be like N~ 2 , 
however this type of function does not fit to over data 
and leads to a big value for the square least error. 

In Figs.^ we present our numerical analysis of the two- 
point GS correlation functions (Eq.(|l4|)) for spin-l-spin-1 
(a,a') = (l,l)andspin-l/2-spin-l/2 (a, a') = (1/2,1/2). 
In Fig.^(a) we observe that without dimerization 7 = 0, 
the spin- 1/2-spin- 1/2 correlation has a slower exponen- 
tial fall-off than the spin-l-spin-1 case which means that 
the quantum fluctuations for spin-1/2 are stronger than 
the spin-1 case. Moreover, they both are positive which 
means that the spins are aligned parallels in each sub- 
lattice (ferro-magnetically) . In Fig.^(b) we make a check 
by comparing with the one-dimensional case. This is 
achieved by using the "snake mechanism" fl3]| : we set 
7=1 and J' = 2 in eq. ([!]) for the SBA configuration. 
Under these conditions, we are dealing with an ordinary 
fcrrimagnetic chain with alternating spins of magnitude 
1/2 and 1 without dimerization and with effective cou- 
pling constant J e // = 2. Thus, as we are plotting the 
connected correlation function 



(14) 



then the case spin 1-1 has opposite sign and happens only 
in one-dimensional limit (J 7 = 2, 7 = 1, SBA). It is due 
to strong quantum fluctuations in the one-dimensional 
case. This is in full agreement with results obtained ear- 
lier in H (see Fig. 5(a) and (b) of this reference.) We 
have reproduced their results in a particular case of our 
study. Moreover, we have extended this numerical anal- 
ysis for real 2-leg mixed spin ladders as explained in the 
other figures. 
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FIG. 7. Correlation functions for both spin (1/2, 1/2) and 
spin (1, 1) (Lanczos): (a) without dimerization; (b)-(d) with 
dimerization (SBA). 

In order to see what is the effect of the interchain cou- 
pling J' on the correlation length £ in the 2-leg ladder as 
compared with the pure one-dimensional case, we plot in 
Fig.^(d) the spin- 1/2-spin- 1/2 case for 7 = and 7 = 1 
at J' = 2 and we observe that in the ladder (7 = 0) there 
is a slower decay than in the chain (7 = 1). This signals 
that the correlation length is bigger in the ferrimagnetic 
ladder than in the corresponding chain: ladder > £,chain- 
Likewise, in Fig. 0(c) we have made a similar analysis to 
see the effect of the dimerization in the staggered config- 
urations (SBA). We observe that in the ladder 7 ^ 1 the 
correlation length £ decreases as the dimerization gets 
stronger. Also, the correlations remain ferromagnetic in 
this spin-1 sub-lattice. 



IV. NUMERICAL PHASE DIAGRAM 

With the numerical tools provided by the Lanczos 
method we can also make an analysis of the phase di- 
agram exhibited by the 2-leg mixed spin ladders in the 
space of couplings of the dimerization 7 and the inter- 
chain coupling constant J'/ J. We have performed this 
analysis for both antifcrromagnetic couplings J' > and 
ferromagnetic couplings J' < (in units of J = 1). More- 
over, we have studied both dimerization patterns, CBA 
and SBA, considered throughout this work. 

In order to set up this phase diagram we determine 
numerically the existence or lack of an energy gap in the 
lowest lying spectrum, as we have done in previous sec- 
tions. For antiferromagnetic J' we find gapless ferrimag- 
netic order for both CBA and SBA configurations. This 
ferrimagnetic order is supported by a non-zero value of 
magnetization. This is in agreement with the approxi- 
mate QRG analysis in ]ll| ] , and thus we do not explicitly 
plot this diagram. The most interesting characteristics 
are found in the ferromagnetic J' region. In this case, we 
find that the GS is a spin singlet while the lowest excita- 
tion is a spin triplet. For the SBA dimerization pattern 
we always find a gapless behaviour for any value of the 
dimerization parameter 7. Since the ground state is a 
singlet the magnetization is zero for all J' < region. 
Thus we conclude that the phase boundary between the 
two gapless phases ( J' < and J' > 0) is at J' = where 
a first-order phase transition occurs. If we consider the 
magnetization as the order parameter then it is discon- 
tinuous at J' = (1-st order transition). The J' = 
transition line modifies the location of phase boundary 
obtained by QRG in the previous study (see Fig. 3 in Ref. 
0)- 

However, the situation turns richer for the CBA dimer- 
ization and this is the case that we plot in FigM 



G 



To obtain this numerical phase diagram we have split the 
range 7 S [0, 1] into 0.1 steps and the range J 1 S [—2, 0) 
into 0.2 steps. Thus, we end up with a 10 x 10 grid of 
numerical points. For each point in this grid, we have 
computed the lowest energy state in the sectors of total 
spin component S z = 0,1. This is done for a series of 
lattice sizes, namely, N = 2 x L with L = 4,6, 8, 10. With 
this data we extrapolate the value of the energy gap for 
infinite lattices and draw a symbols X if the gap is non- 
vanishing and a symbol O otherwise. These extensive 
Lanczos analysis is presented in Fig.^|. 
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FIG. 8. Phase diagram of a 2-leg mixed spin ladder with 
CBA dimerization in the space of couplings 7 v.s. J'/J 
(J' < 0). A cross X symbol denotes a finite gap while an 
O symbol denotes a zero gap. The dashed line is an esti- 
mate for the phase boundary between the gapless and gapful 
phases. 

Here we find two phases depending on whether we are in 
a strong or weak dimerization regime, assuming also that 
the interchain coupling J' is not negligible. We clearly 
detect a gapful phase in the upper part of the diagram 
and a gapless phase in the lower part. From our nu- 
merical data we also estimate the phase separating line 
plotted as a dashed line in Fig.|| and we observe that 
it approaches a horizontal asymptote as the interchain 
coupling becomes very large. As it has been mentioned 
before the whole J 1 < region has zero magnetization 
due to the spin of ground state which is S — 0. So the 
phase boundary between the two gapless phases ( J' < 
and J' > 0) is at J' — accompanying a first-order tran- 
sition similar to SBA configuration. 



V. THE SPIN-PEIERLS INSTABILITY 

We have also carried out a Lanczos study of the spin- 
Peierls (SP) instability. This phenomena has been re- 
alized in several materials and here we anticipate this 
possibility in the ferrimagnetic 2-leg ladders. This tran- 
sition towards a dimerised GS ferrimagnetic ladder is de- 
termined by the competition between the lowering of the 
magnetic energy due to dimerization and the raising of 
the elastic energy due to phonons. According to standard 
terminology, the transition is called unconditional if the 
ground state is dimerised for arbitrary large value of the 
spin-phonon coupling (or small value of dimerization), 
and it is called conditional otherwise. 

Using exact diagonalization techniques it has been pos- 
sible to stablish (l6) that integer spin chains do not have a 
spin-Peierls transition while half-integer spin chains does 
exhibit this transition. This behaviour was predicted 
by Schulz JItJ using bosonization techniques. However, 
other analytical studies had contradicted this conclusion 
Jl6|| . Therefore, the study of the spin-Peierls transition 
provides us with an alternative view of the quantum dif- 
ferences between integer and half-integers spins, in addi- 
tion to the more familiar gapped v.s. gapless behaviour. 

Moreover, it is very interesting to study the SP transi- 
tion for quasi-lD systems as a function of the interchain 
coupling J'. This is usually done within the so called 
chain mean-field theory which underestimates the effect 
of this coupling. Here we can treat its effect without any 
bias. 

AE (y) (1CT 6 ) 





AE„(t) (IP ') 



(c) 
7=0.015 
J'=l 



, AE„(Y) qq ) 
CBA 



7=0.045 




FIG. 9. Plot of the magnetic energy gain (Lanczos) AEo(-j) 
eq. ([[5]) v.s. 1/L for both CBA and SBA configurations 
(spin-Peierls analysis), (a) and (b) for different values of in- 
ter-chain coupling (J') and fixed value of 7 = 0.005. Fixed 
value of J' — 1 is presented in (c) 7 = 0.015 and (d) 7 = 0.045. 
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The magnetic energy gain per bond AEoij) is com- 
puted as follows, 



AE (j) = [E L (0) - E L ( 7 )}/L b , 



(15) 



where L\, = 3L is the number of interacting bonds. It can 
be shown that studying the limit of (|l^) as L — ► oo for 
small 7 reveals the conditional or unconditional character 
of the SP transition. Specifically, this feature is encoded 
in the concavity or convexity character, respectively, of 
the curve of ( p)| ) v.s. 1/L p6[ . Namely, for a finite sys- 
tem the condition for having an unconditional dimerised 
transition is that in the limit of small 7, the energy gain 
( |l5| ) must strictly increase as N increases. This means 
that the stabilization energy always overcomes the elastic 
energy. Otherwise, if the magnetic energy gain does not 
strictly increase with increasing the size, then the transi- 
tion is conditional, i.e., it depends on the stiffness of the 
lattice. This characterization of the conditional v.s. un- 
conditional character of the dimerization transition can 
be recasted in terms of the geometry of the curve for the 
energy magnetic gain (|l^) as a function of the size 1/L: 
when the transition is unconditional the curvature of this 
function is oriented upwards with respect to the x axis 
and we call this a convex curve (see Fig. 1(a)) in p6|); 
while in case of the conditional transition, the curvature 
is oriented downwards with respect to the horizontal axis 
and we call this a concave curve (see Fig. 1(b) in |iq] .) In 
Figs. ^ we present our results. We observe that these 
curves are clearly concave for any values of the J' cou- 
plings that we have investigated. Moreover, this result 
holds true for both CBA and SBA configurations. Thus, 
we find that the SP transition is conditional for a set of 
values in J 1 ranging from to 2.0. We have presented a 
reduced number of plots for the variation of J' in Fig.^ 
because for larger values the curves have a different scale 
and cannot be fitted into one plot. 

We have checked this by taking three different value of 
7 (0.005, 0.015, 0.045). It has been also shown that the 
magnetic energy gain in CBA is bigger than in the SBA 
configuration which is in agreement with Fig.(^j) where 
the CBA configuration is more stable than the SBA one 
against small perturbations. 



VI. CONCLUSIONS 

This paper represents the first extensive numerical 
study using the Lanczos method of quantum ferrimag- 
netism in quasi-one dimensional models with a ladder 
structure. Previous numerical works have dwell upon 
truly one dimensional spin chains with alternating spins 
of magnitude {S 1 ,S 2 ) = (1, 1/2). 

We have presented our results regarding ground state 
properties of a 2-leg ladder with two patterns of dimeriza- 
tion: Columnar (CBA) and Staggered(SBA). Likewise, 



we have seen how the GS properties evolve upon varia- 
tion of the interchain coupling constant J' of the ladder. 

Among the properties we have computed we mention 
the ground state energy, Antiferromagnetic gap, gapless 
behaviour of the Ferromagnetic excitations as we extrap- 
olate our results to the thermodynamic limit, GS corre- 
lation functions and so on and so forth. These results 
have been contrasted with an approximate Spin Wave 
Theory analysis in the linear approximation and found 
a qualitative good agreement in some properties like the 
ground state energy, but not for the gaps. Moreover, 
we have performed a numerical analysis of when a spin- 
Peierls transition can occur and we have reported that 
this transition is conditional. 

We have completed our Lanczos study of 2-leg mixed 
spin ladders with the inclusion of a numerical phase di- 
agram which allows us to clarify the different gapless or 
gapful phases in the space of couplings of the model. 

We believe that these results are of interest for re- 
searchers in the area of quantum spin systems that 
may want to know how the ferrimagnetism of a one- 
dimensional spin chain evolves when an extra chain is 
introduced and treated on equal footing. 
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